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ABSTRACT :  Work  with  finely  ground  AP  led  to  the  infinite  diameter 
relationship 

D^mm/psec)  “  1.146  +  2.576  pQ 

over  the  range  0.55  *  pQ  £  1.0  g/cc.  D  and  pQ  are,  respectively, 
detonation  velocity  and  charge  density.  The  quadratic  mean  error 
found  was  0.040  mm/psec.  Several  computations  of  D^  as  a  function 
of  p0  are  summarized  and  compared  with  the  experimental  data. 
-Estimated  limits  on  the  detonation  temperatures  are  proposed. 

Detonation  velocity  measurements  were  made  on  AP/H.E.  mixtures  con¬ 
taining  of  HMX,  PETN,  and  TNETB.  They  indicated  an  oxidation- 
reduction  reaction  for  AP/HMX  only.  The  AP  in  the  other  two  mixtures 
contributed  to  the  detonation  velocity,  hut  the  additivity  rule  does 
not  seem  applicable  to  mixtures  of  AP  with  organic  H.E. 
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IDEAL  DETONATION  VELOCITY  OP  AMMONIUM  PERCHI ORATE  AND 
ITS  MIXTURES  WITH  H.E. 


In  earlier  work1"^,  very  few  detonation  velocity  (D)  measurements 
were  made  on  ammonium  perchlorate  (AP)  charges  of  loading  density  p  < 

V 

1.0  g/cc  because  such  charges  were  of  very  poor  quality.  Moreover, 
extrapolation  of  the  D  vs  reciprocal  diameter  (d”1)  data  to  the  ideal 
or  infinite  diameter  value  D^^  at  pQ  >  1.0  g/cc  was  suspect  because  of 
the  nature  of  the  D(pQ,  d)  function  for  AP^.  The  present  work  was 
therefore  undertaken  to  obtain  a  D^  vs  pQ  curve  for  pQ  s  l.o  g/cc  and 
to  see  if  it  could  be  extrapolated  to  higher  densities. 

A  new  lot  of  more  finely  ground  AP  was  used  to  prepare  charges  of  . 

0.55  and  0.70  g/cc.  The  results  of  extrapolating  D  vs  d”1  data  were 
then  combined  with  previously  reported  Di  for  0.72  and  1.0  g/cc  to 
obtain  the  desired  D^  vs  pQ  curve  for  pQ  £  1.0  g/cc.  The  extrapola¬ 
tion  of  this  latter  curve  to  higher  densities  was  then  explored  by 
studying  the  D^  values  of  AP/H.E.  mixtures  at  a  high  percent  of  their 
theoretical  density,  i.e.,  in  regions  where  pure  AP  cannot  be  detonated 
at  d  s  7.62  cm. 


EXPERIMENTAL 

A  new  lot  of  AP,  propellant  grade,  was  ground  to  an  average  diameter 
of  11.6  p,  by  surface  mean  (determined  microscopically).  This  ma¬ 
terial,  supplied  by  N.O.S.,  Indian  Head,  Md.,  is  designated.  N-128, 

2i 

According  to  Heywood  ,  the  microscopic  measurement  corresponds  to  a 
Micromerograph  particle  size  measurement  of  0. 67  x  11.6  or  7,8  n.  It 
is  therefore  somewhat  finer  than  the  10-11  p,  (Micromerograph  values) 
size  AP  previously  used.  A  second  lot  of  AP,  prepared  by  N.O.S.  is 
designated  N-133.  It  had  an  averag<  diameter  11. 3  p.  by  surface  mean 
and  a  micromerograph  value  of  7.7  p.  (O.67  x  11.3  =  7,6).  A  third  lot 
of  AP  (N-136)  was  ground  on  the  same  schedule  as  that  used  for  N-133. 
This  lot  had  a  micromerograph  v*lue  of  8.4  p.,  A  fourth  lot  AP  (N-138) 
again  ground  on  the  same  schedule  had  a  micromerograph  value  of  8.9  p.. 
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Charge  preparation  and  handling  were  Identical  to  those  in  previous 
wo;'k^.  AP  charges  of  pQ  s  1.0  g/cc  were  hand  packed  and  fired  im¬ 
mediately  after  preparation.  Charges  were  of  various  diameters,  20.52 
cm  long,  and  boosters  were  50/50  pentolite  (pQ  «  1.56  g/cc)  of  the 
same  diameter  and  5.08  cm  long.  The  experimental  assembly  and  instru¬ 
mentation  (70  mm  smear  camera  with  writing  speed  up  to  4  rmn/p-sec)  were 
also  the  same^. 

Record  reduction  was  carried  out  as  in  the  earlier  work^.  In  addition, 
small  corrections,  not  used  in  the  previous  work,  have  been  made  on 
the  present  data.  They  are  shown  in  the  tables  and  consist  of  the 
following : 

1.  Correction  of  measured  D  to  that  expected  at  pQ  **  1.0  g/cc. 
The  charges  used  in  the  series  of  pQ  ~  1.0  g/cc  had  a  density  range 
of  O.988  to  1.009  g/cc ,  and  the  small  correction  was  made  by  use  of 
the  smoothed  D  vs  pQ  curve  at  the  appropriate  diameter  in  the  earlier 
work.  The  results  in  both  D  and  p0  were  expressed  in  four  figures 
throughout  the  series  of  corrections.  Pinal  results  should  be  rounded 
off  to  three  figures. 

2.  D  was  corrected  for  the  effect  of  using  a  non-planar 
initiating  shock^,  i.e.,  the  measured  velocity  on  the  charge  surface 
was  corrected  to  the  corresponding  axial  velocity.  For  this  correc¬ 
tion,  the  5.08  cm  long  booster  is  treated  as  part  of  the  charge 
length  and  as  if  the  detonation  velocity  of  the  booster  were  the  same 
as  that  of  the  charge.  The  correction  was  made  for  an  axial  point 
3.8  cm  from  the  free  end  of  the  charge ;  this  location  corresponds  to 
the  midpoint  of  that  part  of  the  smear  camera  trace  read  for  the 
velocity  determination. 

3.  This  correction  is  made  for  the  difference  in  booster  and 
charge  detonation  velocities.  It  will  be  fully  described  in  a  sub¬ 
sequent  report.  Its  derivation  assumes  one  constant  velocity  in  the 
booster  and  another  constant  velocity  in  the  charge,  in  other  words, 
no  overboostering .  (For  pentolite,  D  «  7.2  mra/jisec  and  is  always 
greater  than  D  of  AP.) 
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RESULTS  AND  DISCUSSION 
AP 

1516  fine  material,  AP  N-128,  had  a  pour-density  of  about  0.5  g/cc. 
Charges  of  0.55  and  0.70  g/cc  prepared  from  it  were  of  better  uni¬ 
formity  than  the  low  density  charges  in  previous  work,  but  were  still 
of  lower  quality  than  desired.  Because  we  are  now  at  about  the  limit 
of  available  grinding  equipment,  there  seems  little  chance  of  making 
any  further  improvement  in  the  charge  quality  without  using  a  new 
technique  such  as  recrystallization  of  the  AP. 

The  data,  D  as  a  function  of  d  at  various  (0.55  to  1.0  g/cc)  are 

°  l 

shown  in  Table  .1.  Two  sets  of  data  were  pre/iously  reported  ;  two 
are  new.  Each  measured  value  of  D  was  corrected  as  shown  in  the  table. 
Correction  2  (for  non-planar  initiation)  reduces  the  measured  value 
whereas  correction  3  (for  booster  detonation  velocity)  increases  it. 
These  two  corrections  are  therefore  somewhat  compensating  and  the 
corrected  results  will  probably  bracket  the  true  velocity.  Thus 
Table  1  shows  the  minimum  (i >mln,  correction  2  only)  and  the  maximum 
(Dmqy,  corrections  2  and  3).  Each  corrected  D  vs  d  series  was  then 
fit  by  least  squares  to  a  linear  D  vs  d"1  which  could  be  extrapolated 
to  the  infinite  diameter  value  D^.  The  values  so  obtained  and  the 
quadratic  mean  error  (q.m.e.)  for  each  fit  are  also  given  in  Table  1. 
The  average  relative  deviation  for  each  of  the  linear  fits  lies  be¬ 
tween  1.0  and  1.2#. 

Figure  1  displays  the  results  graphically.  The  straight  lines  are 

the  least  square  results  D  vs  d"1,  I'he  data  points  are  the  D_  values. 

max 

(Graphs  based  on  D_.„  look  much  the  same  and  are  so  close  to  the  D 

min  max 

values  that  they  cannot  be  shown  on  the  same  plot.)  The  low  density 
charges  gave  reasonable  linear  D  vs  d"1  curves  despite  the  fact  that 
the  charges  were  not  as  uniform  as  desired.  The  two  curves  for  N-128 
confirmed  a  previous  observation  on  XP-17,  that  the  slope  of  the  D  vs 
d'1  curve  steepens  as  pQ  increases.  In  two  of  the  curves  (XP-11  and, 
to  a  lesser  extent,  N-128  at  pQ  =  0.701  g/cc),  the  value  at  d  =  6.35 
cm  (d”1  «  0.157  cm”1)  seems  to  be  most  out-of-line.  This  is  dis¬ 
cussed  later  in  Appendix  A. 
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TABUS  1  ORIGINAL  AND  CORRECTED  DATA  PC®  SEVERAL  APi 
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TABLE  1  ORIGINAL  AND  CORRECTED  DATA  FOR  SEVERAL  APs.  (Cont'd) 
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1.0  2.0  3.0 

RECIPROCAL  DIAMETER  10  d"1  (cm-1) 


FIG.  1  DETONATION  VELOCITY  VS  RECIPROCAL  DIAMETER  CURVES  FOR  AP 
(Dmax  TABLE  1) 
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The  vs  pQ  data  were  weighted  according  to  the  number  of  measure¬ 
ments  made  at  each  density  and  then  fitted  by  least  squares  to  a 
linear  vs  pQ.  The  results,  given  in  Table  1,  can  be  summarized  as 


Di  (mmAisec)  = 


1.104  +  2.591  P0  min. 
1.146  +  2.576  p  max. 


(la) 

(lb) 


where  pQ  =  0.55  to  1.0  g/cc.  The  curve  for  the  maximum  DA  is  shown 
in  Figure  2.  Extrapolated  results  from  three  different  lots  of  fine 
AP  lie  on  the  same  linear  D  vs  pQ  curve  in  this  experimental  range. 

The  difference  between  the  maximum  and  minimum  values  of  Eqs.  (1)  is 
0.035  and  0.027  ram/iisec  at  the  respective  densities  of  0.55  and  1.0 
g/cc.  It  amounts  therefore  to  1,3  to  0.7#  in  D^.  At  crystal  density 
(1.95  g/cc)  the  difference  is  0.01  mm/p-sec.  Since  the  two  equations 
are  nearly  the  same,  Eq.  (lb)  rather  than  both  will  be  used  for  sub¬ 
sequent  work. 


Eqs.  (l)  are  remarkably  close  to  the  curve  given  by  Andersen  and 
Pesante°  shown  in  Eq.  (2). 

(mm/nsec)  *  1.012  +  2.688  pQ  (2) 

They  differ  markedly  from  the  equations  reported  in  the  earlier 
1  3 

work 


=  -0.016  +  3.784  pQ  0.6sJpo<L.26  g/cc 
and 

=  -0.45  +  4.19  pQ  l«0£po*1.26  g/cc 

although  all  four  equations  agree  at  p_  =  1.0  g/cc.  However,  it  was 

pointed  out  in  our  earlier  work  that  the  derived  value  (pQ  * 

1.26  g/cc)  might  well  be  too  high.  In  view  of  the  slopes  of  Eqs.  (l) 

(with  the  assumption  that  Eqs.  (l)  can  be  extrapolated  above  1.0  g/cc) 

the  first  estimate  of  D.  at  1.26  g/cc  is  too  high  by  about  10#.  It 

1  "51 

was  also  pointed  out  that  the  low  density  charges  were  of  un¬ 
satisfactory  quality  and  hence  only  an  estimate  could  be  made  for 
at  p0  ~  0.6  g/cc.  Again  by  Eqs.  (l),  the  original  estimates  seem 
about  12#  too  low. 
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In  reviewing  the  early  as  well  as  the  present  work  at  low  densities 
we  have  found  some  evidence  that  a  batch  of  AP  will  exhibit  lower 
detonation  velocities,  for  a  given  d  and  pQ,  as  it  ages.  (This  evi¬ 
dence  is  described  in  Appendix  B,  and  additional  data  on  this  topic 
are  being  obtained.)  This  trend  with  age  could  result  from  moisture 
pick-up  by  the  AP*  or  from  particle  agglomeration  or  from  both. 

Since  shots  at  0.6  g/cc  were  made  about  3-6  months  after  the  shots  at 
1.0  g/cc  in  the  earlier  work,  an  aging  effect  of  the  AP  itself  could 
have  contributed  to  the  low  velocities  found. 

An  effort  was  made  to  check  the  present  values  (Eq.  (l))  by  elec¬ 
tronic  measurements  of  D  in  confined  charges.  For  this  purpose  a 
steel  tube  was  used  to  confine  the  charge.  Results  are  shown  in 
Table  2.  Over  the  range  0.55  to  1.0  g/cc,  the  measured  D  is  no  more 
than  2$  lower  than  and  at  1.355  g/cc,  the  difference  is  about  Jjt 
less  than  the  value  from  extrapolation  of  Eq.  (l)  to  1.355  g/cc. 
This  seems  encouraging,  but  at  the  bottom  of  Table  2,  are  given  the 
D  measurements  from  smear  camera  records  of  shots  made  with  unconfined 
charges  of  7.62  cm  diameter.  It  is  evident  that  the  agreen."  .  in  this 
case  is  Just  as  good.  Hence  the  confinement  has  made  the  effective 
diameter  of  the  charge  about  7.62  cm;  it  has  not  extended  our  effec¬ 
tive  experimental  range  in  d  or  solved  our  problems  of  D  vs  d-1 
extrapolation  at  pQ  >  1.0  g/cc. 

D  for  the  shot  at  pQ  =  1.584  (81.2^  TMD)  is  7#  lower  than  the  D^.  ob¬ 
tained  by  extrapolating  Eq.  (lb).  However,  it  is  also  the  highest  D 
value  we  have  observed  in  pure  AP,  and  the  density  is  the  highest  at 
which  we  have  found  pure  AP  detonable.  In  this  case,  the  first  pin 
was  located  1.27  cm  from  the  pentolite  booster.  No  overboostering 
was  observed,  and  the  q.ra.e.  for  this  shot  was  0.213^.  The  sane  AP, 

*Gor’kov  and  Kurbangolina^  have  shown  that  the  critical  diameter  of 
AP  is  much  more  sensitive  to  its  moisture  content  than  is  that  of 
more  common  explosives.  Presumably  moisture  would  also  affect  the 
measured  D, 
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under  the  same  confinement,  fails  at  pQ  =  1.635  g/cc  (83.9#  TMD). 
Nevertheless,  shock  from  the  booster  Induced  some  reaction  which  sub¬ 
sequently  faded.  Velocities  over  the  successive  intervals  between 
pins  were  4.58,  4.62,  4.40,  4.15,  5.50,  2.70,  1.82,  and  1.37  mm/Usec; 
the  10th  pin  was  not  even  shorted  out  by  the  attenuated  front. 

AP  Computations 

In  Table  3  are  the  detonation  parameters  computed  for  a  very  low 
(0.01  g/cc)  density;  the  computed  velocity  is  shown  on  Figure  2. 

O 

Like  other  computed  values  for  very  low  density0,  this  point  fallB 
above  the  extrapolated  ideal  curve  of  Eq.  (lb).  Of  course,  it  has 
been  established  for  some  time  that  current  Ruby  code  computations 
of  the  detonation  velocity  of  AP  do  not  reproduce  the  experimental 
value  at  1.0  g/cc^ 10.  But  the  disagreement  can  be  attributed  to  the 
improper  equation  of  state  or  covolumes  or  both.  Neither  objection 
should  be  applicable  to  the  present  computation  made  at  such  a  low 
density  that  ideal  gas  behavior  can  be  assumed.  Moreover,  this  cal¬ 
culation  gives  us  some  information  about  the  detonation  reaction  of  AP. 

The  low  pressure  burning  of  A?  has  been  extensively  studied  and  it  is 
known  that  N20  or  NO  (depending  on  the  temperature)  is  a  major  product. 
However,  in  detonation  at  low  pressure  (low  charge  density)  we  know 
from  the  data  of  Table  3  that  the  decomposition  reaction  approximates. 

NH4C104  — 0.5  N2  +  HC1  +  1.5  H20  +  1.25  0g  (3) 

and  that  only  a  small  amount  of  Cl  is  formed  at  0.1  kbar  and  1957°K, 
which  is  probably  near  the  maximum  temperature  to  be  expected  at  any 
pressure.  Eq.  (3)  is,  in  fact,  the  limit  equation  for  "high  pressure, 
high  temperature"  reaction  which  Simchen'1'1  derived  from  thermal  de¬ 
composition  data.  In  the  practical  range  of  densities,  say  from  0.55 
to  1.0  g/cc,  the  reaction  pressures  in  detonation  will  be  much  higher, 
however.  If  we  use  the  data  of  Table  1  and  the  approximation 
Pj  ^  0.25  p0  the  range  is  9-54  kbar,  i.e.,  a  100  fold  or  greater 

increase  in  pressure.  Under  these  circumstances,  a  shift  in  the  de¬ 
composition  reaction  toward 

nh4cio4 — -0.5  n2.  +  0.5  Cl2  +  2  h2o  +  o2  (4) 
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TABLE  3 

COMPUTED*  VALUES  FOR  A?  AT  p0  *  O.Ol  g/cc 


C-J  Values 

Velocity 

1.525  nun/p-sec 

Pressure 

104.6  bars 

Temperature 

1957°K  (chemical  energy  342  cal/g) 

Density 

0.0180  g/cc 

Gamma 

1.245 

Total  gas 

36.02  m.oles/1000g  AP 

Detonation 

Moles/ 

Products* ** 

lOOOg  AP 

Cl 

0.269 

C!g 

0.438 

H20 

15.558 

HC1 

7.566 

Ng 

4.201 

NO 

0.107 

°2 

10.297 

*  Ruby  code  with  ideal  gas  law  as  equation  of  state 
for  the  products. 

**  In  decreasing  amounts:  NOg,  Hg,  NgO,  and  NH,  with 
jcaximum  concentration  (NOg)  of  0.0022  moles/iOOOg  AF. 
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would  be  favored.  It  will  be  ahown  a  little  later  that  equilibrium 
computations  do  indicate  that  Eq.  (4)  describes  the  detonation  reac¬ 
tion.  It  is  also  worth  noting  that  the  first  mechanism  releases 
slightly  less  energy  (346  vs  405  cal/g)  and  slightly  more  gas  (36  vs 
34  moles Ag)  than  the  second. 

It  was  remarked  above  that  the  current  Ruby  code  does  not  produce  the 
correct  value  for  AP.  This  is  true,  but  it  does  not  seem  to  be 
exclusively  the  fault  of  the  equation  of  state  used,  and  the  current 
computations  do  seem  to  predict  the  correct  relative  trend  of  Di  vs 
pQ,  Consequently,  it  seems  worthwhile  to  review  briefly  those  com¬ 
putations  which  have  been  made  on  AP  and  determine  what  information 
they  give. 

The  Ruby  code  uses  the  BKW  equation  of  state 

PV  /RT  -  1  +  xepx  where 

o 

x  -  xEVjkj/7  ( T  +  0)a 

and  P,  Vg,  R,  T,  v^k3  are  pressure,  molar  volume  of  gas  products, 
gas  constant,  absolute  temperature,  mole  fraction  of  gas  product  i, 
and  covolume  of  gas  product  i,  respectively.  The  “constants"  or 
equation  of  state  parameters  to  be  assigned  values  are  the  set  (x, 
a,  0,  0)  and  the  covolumes  k^.  This  code  solves  the  hydrodynamic  and 
the  thermodynamic  relationships  simultaneously  by  an  iterative  method. 
Its  computed  v.^  are  therefore  correct  amounts  of  equilibrium  products 
for  the  conditions,  Pj,  Tj,  it  has  computed  from  the  hydrodynamic 
conditions;  of  course,  P^  and  T^  can  be  wrong  to  the  extent  that  the 
equation  of  state  is  inadequate. 

Table  4  contains  the  results  of  the  computations  for  the  parameters 
(II.85,  0.5,  0.09,  400)9  and  k±  set  I9  and  (12.685,  0.50,  0. 09585, 
400)12,‘1'^  and  k^,  set  II12* Although  the  parameters  have  been 
changed  slightly  and  the  set  of  k^,  appreciably,  the  results  of 
Table  4  show  little  difference  between  the  values  produced  by  these 
two  calculations.  In  particular,  both  show  vi  closely  approximating 
Eq.  (4)  over  the  range  of  P^  of  17  to  187  kbar  and  of  780°-1722°K. 
Moreover,  both  overestimate  D^,  (and  hence  Pj)  although  they  exhibit 
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the  expected  trend  of  D,  increasirg  with  prt.  In  particular,  our  Ruby 

IP  1  ° 

computations  produce  about  the  same  slope  as  that  found  for  the 
experimental  curve  (See  Figure  3)  although,  the  absolute  value  of 
Is  about  0.7  mm/psec  too  high.  In  view  of  the  correct  slope,  it  is 
possible  that  the  relative  trend  is  fairly  good. 

A  much  earlier  computation  was  made  with  the  BKW  equation,  the  set  of 
parameters  (1.0,  0.25,  0.30,  0) 1/1,15  and  the  k ±  set  III.14,15  In  this 
case,  the  equilibrium  was  not  calculated;  instead  Eq.  (4)  was  assumed. 
(Our  current  computations  seem  to  Justify  this  assumption  completely; 
hence  any  difference  in  results  must  be  attributed  tc  the  combination: 
different  parameters  and  different  k^.)  Table  5  and  Figure  3  show 
the  results  of  this  earlier  computation,  which  approximate  the  experi¬ 
mental  values  much  better  than  the  current  calculations,  despite  the 
fact  that  the  BKW  equation  is  still  used. 

Finally,  another  of  the  earlier  computations,  which  used  the  arbitrary 
decomposition  mechanism*  of  Eq.  (4)  was  carried  out  with  an  Abel 
equation  of  stated  This  too  is  shown  in  Figure  3  and  is,  in  fact, 
the  closest  to  the  experimental  curve.  However,  we  attribute  this 
closest  approach  not  to  a  superior  equation  of  state,  but  to  the  fact 
that  these  computations  were  adjusted  to  agree  with  the  experimental 
point  at  1.0  g/cc10.  A  similar  adjustment  to  the  Reference  14  cal¬ 
culation  might  well  have  yielded  Just  as  good  agreement  with  the 
experimental 

It  is  not  surprising  that  parameters  of  almost  any  real  gas  equation 
of  state  can  be  adjusted  to  give  agreement  for  any  particular  ex¬ 
plosive,  but  it  is  quite  interesting  that  the  parameters  of  Reference  14 
were  not  so  adjusted.  They  were  the  same  parameters  which  had  up  to 
that  time  given  acceptable  agreement  for  conventional  organic  ex¬ 
plosives,  This  raises  the  suggestion  that  the  refinement  of  the  com¬ 
putation  in  the  present  Ruby  codes  might  even  be  a  step  backwards. 

Although  with  the  proper  choice  of  parameters  and  covolumes  either 
the  BKW  or  the  Abel  equation  can  be  used  to  predict  (and  pre¬ 
sumably  Pj)  in  agreement  with  those  obtained  experimentally, 

wThe  equilibrium  between  Cl  and  Cl2  was  considered,  but  as  we  have 
already  seen,  this  has  little  effect  on  the  results. 
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TABLE  5 

COMPUTED*  DETONATION  PARAMETERS  FOR  AP 
Reference  (14) 


oj 

ooj 

Q.'vJ 

D 

mm/usec 

P 

kbar 

JL 

0.8 

3.35 

25 

1940 

1.0 

3.82 

39 

1880 

1.2 

4.51 

57 

1800 

1.4 

4.80 

79 

1710 

1.6 

5.32 

106 

1620 

1.95 

6.29 

118** 

1400 

*  Covoluraes  used  were  those  of  Reference  15  and 
514,  277  for  Cl2,  02. 

**  This  is  obviously  a  typographical  error  and  all 
of  the  values  computed  at  this  density  have  been 
discarded. 
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(1)  EXPERIMENTAL,  EQN.  (lb) 

(6)  REF.  (6),  COOK  EQN.  ADJUSTED  TO  AN 
EXPERIMENTAL  VALUE,  ARBITRARY 
DECOMPOSITION. 

(9)  REF.  (9),  BKWEQN.  (11.85,0.5,0.09,400) 

(12)  REF.  (12),  TWO  POINT  COMPUTATION 
BKW  EQN.  (12.685,0.50,0.09585,400) 

(14)  REF.  (14),  BKWEQN.  (1.0,0.25,0.30,0) 
ARBITRARY  DECOMPOSITION 


LOADING  DENSITY  p0  (g/cc) 


FIG.  3  COMPARISON  OF  EXPERIMENTAL  Dj  VS  p0  CURVE  FOR  AP 
WITH  SEVERAL  COMPUTED  CURVES 
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temperature  is  another  matter.  The  fundamental  objection  to  the  Abel 
and  the  ideal  gas  equation  in  any  practical  range  of  detonation  of 
solid  H.E.  is  that  they  include  no  intermole cular  forces;  such  forces 
must  be  operative  in  real,  dense  gases.  Hence  all  the  energy  pro¬ 
duced  by  the  detonation  appears  as  increased  temperature  of  the  gas 
products,  and  the  calculated  T^  will  be  too  high.  The  fundamental 
objection  to  the  BKW  equation  of  state  is  that  it  assumes  only  re¬ 
pulsive  intermole cular  forces  and  hence  much  of  the  energy  appears  in 
a  density  dependent  term  of  the  internal  energy  of  the  products  and 
does  not  contribute  to  raising  their  temperature.  Hence  this  equation 
of  state  should  predict  T^  which  are  too  low. 

It  seems  reasonable  to  take  about  2000°K  as  an  upper  limit  of  T^  for 
AP  at  any  charge  density.  This  is  based  on  the  computed  Tj  of  1932°K 
for  low  pressures  where  Eq.  (3)  is  operative,  allows  for  some  increase 
as  the  decomposition  shifts  to  Eq.  (4)  and  its  15$  increase  in  chemi¬ 
cal  energy,  but  also  keeps  the  net  temperature  increase  low  because 
the  repulsive  forces  also  become  more  important  with  increasing 
pressure  and  density.  The  trend  of  decreasing  T^  with  increasing  pQ 
can  be  accepted  as  a  reflection  of  increasing  molecular  repulsion. 

The  adiabatic  flame  temperature  (T^)  at  atmospheric  pressure  is 

*  1 
chosen  as  a  lower  limit  for  possible  Tj . '  Prom  previous  computations, 

generally  Tf  <  T^  (95$  TMD)  <  Tj  ($  TMD  <  95)  when  the  Tj  were  com¬ 
puted  with  the  BKW  equation  of  state.  Prom  the  arguments  Just  given 
above,  these  Tj  values  are  already  too  low.  Hence  Tf  seems  a  very 
conservative  lower  limit  for  the  true  detonation  temperature.  For 
AP,  Tf  is  1375°K.17  Moreoyer,  1375°K  is  only  about  100°K  above  the 
minimum  temperature  required  for  the  lower  pressure  limit  of  steady 
state  burning  of  AP*.  As  the  pressure  is  increased,  Tf  also  in¬ 
creases  . 

It  is  pertiaps  a  matter  of  interest  that  Belyaev  and  Lukashenya  have 
reporteathat  explosives  burning  at  one  atm.  exhibit  flame 

*E.  A.  Arden,  J.  Powling,  and  W.  A.  W.  Smith,  "Observations  on  the 
Burning  Rate  of  Ammonium  Perchlorate",  Combustion  and  Flame,  6,  21 
(1962). 
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temperatures  below  the  computed  values.  Their  measurements  were  made 
with  thermocouples  and  checked  with  computations  based  on  the  tem¬ 
perature  coefficient  of  the  burning  rate.  However,  Pokhil  and 
Mal'tsev1^  used  an  optical  method  to  measure  flame  temperatures  in 
the  range  20-100  atm.  Their  results  at  20  atm.  on  tetryl,  PETN,  and 
RDX  were  very  close  to  our  computed  flame  temperatures1^  (p  «  1  atm.). 

Figure  4  shows  these  two  limits  for  the  detonation  temperature  and 
the  variation  Tj  vs  pQ  from  the  various  computations  already  compared 
for  Dj  in  Figure  3.  Only  one  of  the  four  curves,  that  of 
Reference  (14)  stays  within  the  two  limits  and  that  would  have  to 
exhibit  a  small  reversal  to  extrapolate  to  the  low  pressure  value. 
Reference  (9)  and  Reference  (12)  values  are  much  too  low,  and 
Reference  (6)  values  seem  too  high  although  that  curve  too  might, 
with  a  small  reversal,  extrapolate  to  the  low  pressure  value. 

AP/HE 

Investigation  of  the  model  AP/HE  is  a  portion  of  our  planned  systematic 
study  of  propellant  models.  The  high  density,  high  #TMD,  part  of  the 
experimental  range  was  investigated  at  this  time  with  the  hope  that 
AP/HE  chargee,  d  < 7.6  cm  and  #  TMD  >  95#,  could  be  initiated  to 
detonation,  and  that  the  validity  of  extrapolating  Eq.  (lb)  to  pQ  >  1 
g/cc  might  be  tested  by  use  of  the  additivity  rule.  The  first  hope, 
i.e.,  detonation  at  high  #  TMD,  was  realized,  but  the  second  was  not, 
as  the  data  below  will  show. 

A  number  of  preliminary  shots  were  made  and  these  data  are  reported 
in  Appendix  C.  We  found  that  the  90/10  mixtures  seemed  less  uniform 
than  the  80/20;  the  finer  the  AP,  the  higher  the  measured  D;  the 
finer  the  H.E.,  the  more  detonable  the  charge.  As  a  result  of  this 
exploratory  work,  we  decided  to  work  with  80/20  mixtures  of  as  fine 
AP  and  as  fine  H.E.  as  were  available.  For  all  of  the  AP/H.E.  work, 
charges  were  22.9  cm  long  and  the  1.56  g/cc  pentolite  boosters  (of 
the  same  diameter  as  the  charges)  were  2.5  cm  long. 

To  assess  the  extent  of  the  effect  on  the  detonation  velocity  of  the 
reaction  of  the  component  AP  in  AP/H.E.  mixtures,  it  is  desirable  to 
have  data  on  analogous  Inert/H.E.  mixes.  However,  there  are  very  few 
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data  for  mixtures  containing  80#  inert.  Reference  (20)  contains  data 
for  PETN/Selectron*,  20/80;  and  reference  (21),  data  for  15  p.-PETN/ 
polynitroure thane**.  At  the  common  diameter  of  41.35mm,  both  refer¬ 
ences  give  D  «  5.5  .vUsec.  Hence  we  know  that  these  20/80  mixtures 
made  with  fine  PETN  can  detonate.  Moreover,  reference  (20)  shows 
that  inert  diluents  lower  the  detonation  velocity  according  to  their 
density  and  the  volume  fraction  that  can  be  incorporated  into  the 
mixture;  the  higher  the  density  or  the  volume  fraction,  the  greater 
the  amount  by  which  D  is  reduced.  Finally,  reference  (21)  shows  that 
the  diameter  effect  on  D  for  20/80,  PETN/Inert  is  relatively  large. 

To  obtain  a  more  quantitative  estimate  of  the  effect  of  an  inert,  a 
few  5.08  cm  diameter  charges  of  HMX/Wax  and  HMX/KCt,  20/80,  were 
tested.  The  data,  given  in  Table  6,  show  that  HMX/Wax  behaves  as  a 
group  1  explosive  and  has  a  relatively  high  critical  value  of  the  % 

TMD  below  which  the  charge  cannot  be  detonated.  They  also  show  that 
KCt,  the  higher  density  inert,  lowers  the  detonation  velocity  more 
than  wax  although  the  volume  fractions  are  0.86  (wax)  and  0.77  (KCt). 
XCt  (1.99  g/cc)  is  a  better  match  for  inert  AP  (1.95  g/cc)  than  wax 
(1.0  g/cc). 

In  order  to  obtain  an  estimate  for  comparable  PETN/Inert  mixes,  we 
assume  that,  at  any  fixed  charged  diameter  and  fixed  weight  fraction*** **** 
of  the  same  diluent,  the  ratio  of  the  detonation  velocities  of  the 
analogous  PETN  and  HMX  mixtures  will  be  the  same  as  the  ratio  of  the 
D±s  of  the  undiluted  H.E.,  i.e.,  (8.456/9.140)  -  0.925.*#**he  esti¬ 
mates  for  20/80  PETN/Inert  at  d  *  4.08  cm  are  4.74  and  4.57  mm/iisec 
for  wax  and  KC-t  respectively.  These  values  and  those  of  Table  6  for 
the  HMX  mixtures  can  be  compared  directly  to  the  measured  Ds  of  5.08  cm 

*  A  thermosetting  resin  of  1.20  g/cc  density. 

**  Density  of  1.21  g/cc. 

***  For  PETN  and  HMX,  the  same  weight  fraction  of  the  same  inert 

diluent  results  in  nearly  the  same  volume  fraction.  The  small 

20 

amount  of  porosity  will  also  have  some  effect  ,  but  it  has  been 
ignored  in  these  estimates. 

****  Reference  22  curves  used  for  PETN  at  p0  =  1.78  g/cc  and  for  RDX 
at  p0  «  1.90  g/cc. 
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TABLE  6 


DATA  FOR  HMX/INERT,  20/80 
(d  «  5.08  cm) 


Shot  Ho, 


g?8c  g+MD 
HHX/Waxa 


D 

gp/useo 


516  0.543h  49.1 

517  0.557h  50.4 

519  0.697h  63.1 

520  0.901H  83.4 

521  1.0801  97.7 

fflg/*Clb 


F 

F 

5.12 


fi  >  8.8  cm 
6  >  15*3  cm 


1.9191  97.5  4.716 

1.9201  97.5  4.715 


a.  pv  «  1.105  g/cc.  Voidless  densities  of  K.4X  and 
Wax  are  I.90  and  1.00  g/cc,  respectively. 

b.  pv  =  1.971  g/cc.  Voidless  density  of  KC1  is  1.99 
g/cc. 

HMX,  X586,  is  Grade  II,  Class  E  with  an  average  particle 
size  of  about  1%. 
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diameter  charges  of  H.E./AP,  20/80.  However,  if  an  estimate  of  D.^ 
rather  than  D(d  »=  5.1  cm)  is  desired,  we  can  use  the  factor  of 
[Di/0(d  -  5.1)]  »  1.06  from  reference  (21).  Estimates  of  Dj^  for  the 
H.E./KC4,  20/80,  mixes  are*  then  5.0  and  4.6  mm/psec  for  HMX  and  PETN, 
respectively. 

HMX  and  PETN  mixtures.  Both  the  fine  HMX  and  PETN  are  stored  uncler 
water  and  must  be  dried  before  they  are  used  to  prepare  charges.  No 
method  of  vacuum  drying  was  entirely  successful  in  preventing 
agglomeration.  Consequently  the  charges  were  less  uniform  than 
desired  and,  in  some  series,  twice  the  usual  number  of  shots  were 
made  to  compensate  somewhat  for  the  variation  in  charge  quality. 

The  data  for  AP/HMX,  80/20,  at  95.4*  TMD  are  given  in  Table  7  ano 
plotted  in  the  top  curve  of  Figure  5*  The  similarity  of  this  D  vs 
d  curve  to  those  of  AP/Wax,  J  in  which  an  oxidation-reduction 
reaction  occurs,  shows  that  the  HMX  is  acting  as  a  fuel  as  well  as  an 
explosive.  HMX  is  C^HgNgOg  and  has  an  oxygen  balance  of  0*  to  CO  and 
of  -22*  to  C02.  The  curve  of  Figure  5  indicates  that  the  effect  of 
the  fuel-oxidant  reaction  on  D  is  not  evident  until  larger  diameters 
(say  above  5  cm)  are  reached.  This  was  also  evident  in  the  waxed  APs 
and  suggests  that  the  oxidation  reduction  reaction  requires  a  longer 
time  than  decomposition  of  H.E.  Inasmuch  as  this  reaction  rate  must 
be,  to  some  extent,  diffusion  controlled,  the  longer  time  and  reaction 
zone  are  to  be  expected. 

If  the  2-polnt  extrapolation  of  results  for  the  two  largest  charges 
of  Figure  5  is  made,  the  resultant  value  is  7.49  mm/usec.  For 
the  assumptions  of  no  oxidr,cion-reduction  reaction  and  simple  addi¬ 
tivity  of  the  detonation  rates  of  AP  and  HMX,  the  computed  D^  of  AP 
from  the  measured  value  for  the  mixture  is  7.15  mm/iisec.  The  Di 
value  for  AP  at  the  same  *  TMD  (95.4*)  by  extrapolation  of  Eq.  (lb) 
is  only  5.94  mm/psec.  This  is  additional  evidence  that  an  oxidation 
reaction  did  occur  and  contributed  to  the  measured  D  of  the  larger 
charges  of  AP/HMX,  80/20,  Finally,  the  three  AP/HMX  series,  reported 
in  Appendix  C,  also  show  excessively  high  D  values.  We  conclude  that 
the  detonation  products  of  HMX  are  oxidized  by  those  of  AP,  and  that 
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TABLE  7 

DATA  FOR  AP/HMX,  80/20 
AT  95. 4g  TMD 

Charge 


Shot  No. 

Dlam. 

d,cm 

S%c 

D 

mm/p.sec 

Dl 

D2 

525 

2.540 

1.8491 

6.225 

6.214 

6.215 

524 

3.495 

1.8491 

6.249 

6.229 

6.229 

526 

5.080 

1.8541 

6.486 

6.442 

6.443 

533 

7.620 

Av. 

1.8511 

1.851 

6.896 

6.791 

6.793 

AP,  N133  (7.7m-  by  micromerograph) 

HMX,  X586  (1%) 

Voidless  density  of  AP/HMX,  80/20,  Is  1.940  g/cc 
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this  contribution  to  the  detonation  behavior  will  also  occur  in 
propellants  containing  HMX. 

In  contrast  to  HMX,  PETN  (O^HgN^O^g)  has  an  oxygen  balance  of  -10#  to 
C02  and  should  be  a  much  poorer  fuel.  Data  for  AP/PETN,  80/20  appear 
in  Table  8,  and  are  plotted  in  the  lower  curve  of  Figure  5.  Here 
there  is  no  indication  of  oxidation-reduction  reaction,  and  the  least 
squares  fit  of  the  D  vs  d”1  data  extrapolates  to  Di  =  6.09  mra/psec  as 
compared  to  the  estimated  D^  of  4.6  mm/psec  for  the  analogous  AP/KC-t 
mixture.  Hence  there  seems  no  doubt  that  the  AP  is  contributing  to 
the  detonation  velocity  of  the  AP/PETN  mixture  of  99.1#  TMD  although 
AP  itself  will  not  detonate  at  that  low  porosity.  From  the  D.^  of  the 
PETN  mixture,  we  can  compute  by  simple  additivity  D^  for  AP  at  99.1# 
TMD  (1.931  g/cc)  as  5.51  mm/psec.  Eq.  (lb),  on  the  other  hand,  gives 
Di  (Po  =  s/cc)  as  6.12  mm/pLsec  or  10#  higher.  There  are  two 

possibilities:  (1)  that  the  AP  in  a  mixture  with  common  organic  H.E. 
cannot  contribute  fully  to  the  detonation  velocity  of  the  mixture  at 
the  relatively  small  charge  diameters  (d  <  7.6  cm)  used  in  the  present 
work  and  (2)  that  the  Dj  vs  pQ  curve,  Eq.  (lb),  cannct  be  extrapolated 
to  pQ  >  1  g/cc  because  the  ideal  curve  for  AP  falls  below  this 
straight  line  at  the  higher  #  TMD.  No  selection  can  be  made  between 
the  alternate  possibilities  without  working  at  larger  charge  diameters 
than  our  firing  facilities  permit.  Our  opinion  is  that  the  first 
possibility  is  the  more  probable,  and  that  the  additivity  treatment 
is  not  applicable  to  AP/H.E.  mixtures  (at  least  in  relatively  small 
diameter  charges)  because  the  reaction  time  of  the  AP  in  detonation 
is  very  much  larger  than  that  of  the  organic  H.E. 

It  is  of  interest,  though  possibly  fortuitious,  that  a  2-point  ex¬ 
trapolation  of  the  lower  and  presumed  non-reactive  portion  of  the 
AP/HMX  curve  (Figure  5)  gives  a  D^  value  of  the  mixture  (95.4#  TMD) 

as  6.26  mm/p-sec.  By  additivity,  this  leads  to  D.^  (1.86  g/cc)  of 

5.51  mm/p sec  for  AP  as  compared  to  the  extrapolated  value  of  5.94 
mm/psec  from  Eq.  (lb).  In  other  words,  the  AP  in  AP/HMX  charges  of 

d  <  4  cm  and  in  AP/PETN  at  d  <  7.6  cm  seems  to  be  making  the  same 

contribution  to  the  D  of  the  mixture. 
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TABLE  8 

DATA  FOR  AP/PETN,  80/20 
AT  99.1#  TMD 


Shot  No. 

d 

cm 

8° 

g/cc 

D 

Di 

Dp 

D2* 

mm/iisec 

Po*1.895 

528 

2.540 

1.895 

5.771 

5.761 

5.811 

5.76I 

5.761 

527 

3.495 

1.890 

5.850 

5.812 

5.826 

529 

5.080 

1.898 

5.992 

5.951 

5.953 

5.944 

532 

6.350 

1.877 

5.901 

5.858 

5.845 

5.894 

565 

2.540 

1.906 

5.677 

5.667 

5.668 

5.657 

564 

3.495 

1.909 

5.748 

5.730 

5.691 

565 

5.080 

1.908 

5.917 

5.879 

5.842 

551 

7.370 

1.873 

6.010 

5.924 

5.931 

5.993 

Di  (Po  = 

1.895)  * 

6.09  (1.1#)  -  1.04  (25#) 

d-1 

(0.068)  1 

*  Small  corrections  for  density  made  with  slope  of  2.82 
mu  iisec**'1'  cc  g"1.  This  was  derived  from  slopes  of  the 
Di  vs  Po  curves  for"  Ap»  Eq.  (lb),  and  PETN22;  the 
respective  slopes  used  were  2.576  and  5.880,  The  de¬ 
rivation  assumes  simple  additivity  effect  on  the  detona 
tion  velocity  of  the  mixture. 


AP,  N153 

PETN,  X591>  Primer  grade,  supplied  by  Hercules  Powder  Co. 
Voidless  density  of  AP/PETN,  80/20,  is  1.915  g/cc 
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TNETB  and  AP/TNETB.  Although  no  evidence  of  an  oxidation- reduct Ion 
reaction  was  found  in  the  AP/PETN  mixtures,  we  wanted  to  check  that 
result  with  an  H.E.  containing  more  oxygen.  The  only  such  explosive 
available  in  sufficient  amount  was  TNETB,  CgHgNgO^,  with  an  oxygen 
balance  to  C02  of  -4#.  Our  supply  of  TNETB  was  stored  dry  in  the 
magazine,  and  there  was  enough  in  moderately  fine  granular  form  to 
prepare  H.E.  charges.  But  for  preparing  AP/TNETB  mixtures  some  of 
the  explosive  had  to  be  recrystallized  to  obtain  a  fine  material  of 
about  10  p  particle  size. 

In  order  to  interpret  velocity  data  obtained  on  mixtures  of  AP/TNETB, 
it  was  necessary  to  establish  the  vs  pQ  relation  for  TNETB  Itself. 
The  equation  currently  used  is  that  reported  in  NOLTR  65-218  and  is  in 
error.*  Fortunately,  the  error  appears  to  be  small. 

Table  9  shows  the  shots  fired  in  the  present  work  and  derivation  of 
the  value  (pQ  *  1.76  g/cc)  **  8.389  mm/psec  for  the  isostatically 
pressed  charges.  The  data  and  the  extrapolation  are  illustrated  in 
Figure  6  which  emphasizes  the  very  small  diameter  effect  on  D  for  this 
material  at  high  density  (e.g.,  the  value  at  d  *  2.5  cm  is  within  1# 
of  D^) .  This  Justifies  the  approximations  of  Table  10  where  D 
measured  at  diameters  of  5  and  4  cm  for  densities  of  I.69  and  1.29 
g/cc,  respectively,  are  considered  good  approximations  to  the  cor¬ 
responding  D^.  In  fact,  the  subsequent  treatment  Indicates  that  D 
measured  at  d  =  1.28  cm  and  pQ  =  1.27  g/cc  is  also  a  good  approximation 
to  Di# 

Before  leaving  Figure  6,  it  should  be  noted  that  two  D  values  from 

hydraulically  pressed  charges  (p  ~  1.76  g/cc)  lie  appreciably  below 
-T  0 

the  D  vs  d  curve  derived  for  data  obtained  from  isostatically 
pressed  charges.  The  same  difference,  i.e.,  lower  velocities  measured 
on  charges  prepared  in  the  hydraulic  press  than  on  those  of  the  same 
density  from  the  isostatic  press,  has  also  been  noted  and  reported  in 
recent  work  on  AP  and  nitroguanidine ,  Not  enough  work  has  been  done 
here  or  elsewhere  to  show  whether  the  same  extrapolated  value  (Di) 
could  be  obtained  with  the  two  methods  of  preparation.  Charges  for 
D  vs  d”1  curve-j  should  all  be  prepared  the  same  way. 

*  N.  L.  Coleburn,  private  communication. 
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Table  10  shows  the  data  which  were  assembled  to  derive  the  D,  vs  o 

i  ro 

curve.  They  give 

D^  (mm/psec)  =  1.947  +  3.660  pQ,  o  =  0.015  (5) 

where  pQ  is  in  g/cc.  This  is  remarkably  close  to  the  most  recent 
curve  reported  for  PETN*  and  not  very  different  from  the  TNETB  curve 
reported  in  NOLTR  65-218.**  At  high  densities  the  latter  curve  and 
Eq.  5  produce  essentially  the  same  results. 

The  curve  for  Eq.  5  and  the  data  points  are  shown  in  Figure  7  (solid 
line)  together  with  three  computations  of  D1  vs  pQ  (dashed  lines)  for 
TNETB.  The  Ruby  code  result  with  the  RDX  parameters  gives  a  better 
approximation  than  the  code  with  its  original  parameters,  but  the 
best  approximation  is  that  obtained  with  the  parameters  of  Snay  and 
Christian  (See  NAVORD  2611). 

The  AP/TNETB,  80/20,  charges  appeared  more  uniform  than  those  pre¬ 
pared  from  either  of  the  other  two  mixtures.  (A  barrel  blender  was 
used  to  prepare  the  entire  batch  of  AP/TNETB  in  one  operation.  This 
may  have  improved  the  uniformity  of  this  mixture  over  that  attainable 
in  the  piece-wise  mixing  used  for  earlier  batches.)  Nevertheless 
the  charge  surfaces  were  speckled  with  agglomerates  of  up  to  0.8  mm 
diameter.  The  data  obtained  on  a  diameter  series  for  this  mixture 
appear  in  Table  11.  With  the  exception  of  the  5  cm  dlam  shot,  for 
v/hlch  the  value  seems  too  high,  all  D's  in  this  series  fell  within 
the  range  of  6.01  to  6.09  mm/p-sec.  Apparently  the  diameter  effect  on 
D  is  so  small  for  this  mixture  (possibly  as  small  as  that  for  TNETB 
alone  (Figure  6)  that  it  is  not  detectable  above  the  small  random 
errors  introduced  by  non- uniformity  of  charge,  small  errors  in  density, 
and  the  usual  record  reading  errors.  Consequently  the  values  of 
Table  11  have  been  averaged  to  obtain  a  of  6.07  mm/p.3ec  for  AP/ 
TNETB,  80/20,  at  97.6$  TMD. 

*  Di  =  2.00  +  3.58  pQ,  E.  L.  Lee  and  H.  C.  Hornig,  "Equation  of  State 
of  Detonation  Product  Gases,"  in  Preprints  of  12th  Symposium  (Inter¬ 
national)  on  Combustion,  (1968). 

**  B±  =  1.597  +  3.850  pQ 
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.  used  for  unconfined,  pressed  charges  only.  Additional  measurements 
cast  charges  and  on  very  porous  charges  supported  in  glass  or  plastic 
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LOADING  DENSITY  pQ  (g/cc) 

FIG.  7  D;  VS  pQ  FOR  TNETB:  EXPERIMENTAL  AND  COMPUTED  CURVES 


vs  d  DATA  FOR  AP/TNETB,  80/20,  AT  97. TMD 
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The  AP  in  the  mixture  appears  by  the  additivity  rule  to  have  a  of 
5.5  mm/psec  at  pQ  =  1.904.  Eq.  (lb),  on  the  other  hand,  gives  a 
value  of  6.05  mm/psec  or  S#  higher.  Hence  the  AP  in  this  mixture  is 
behaving  in  much  the  same  way  it  did  in  the  analogous  PETN  mixture; 
i.e.,  it  reacted  sufficiently  to  make  a  contribution  to  the  detonation 
velocity  of  the  mixture,  but  the  contribution  was  not  as  large  as 
that  indicated  by  extrapolation  of  Eq.  (lb). 

SUMMARY 

Finely  ground  APs  were  used  to  derive  the  infinite  diameter  detonation 
velocity  -  loading  density  relation 

D.^  (mm/psec)  =  1.146  +  2.576  pQ  (lb) 

over  the  range  0.55  <  pQ  <  1.0  g/cc.  The  highest  velocity  actually 
measured  was  4.86  mm/psec  observed  in  a  confined  charge  at  pQ  =  1.584 
g/cc  (81.2$  TMD) .  Various  hydrodynamic-thermodynamic  computions  re¬ 
produce  the  slope  of  Eq.  (lb),  but  some  give  values  that  are  too 
high  by  0.7  mm/psec.  None  of  these  computations  produce  reasonable 
values  of  the  detonation  temperature  which  is  estimated  to  lie  be¬ 
tween  2000°K  and  1375°K,  the  adiabatic  flame  temperature. 

Mixtures  of  AP/H.E.,  80/20,  detonate  at  99$  TMD.  Results  from  diameter 
series  give  values  for  mixtures  with  HMX  indicating  an  oxidation- 
reduction  reaction.  Although  neither  PETN  nor  TNETB  mixtures  appear 
to  have  such  a  reaction,  the  contribution  of  AP  to  their  D^  is  less 
than  that  Indicated  by  Eq.  (lb).  It  is  suggested  that  the  additivity 
rule  is  not  applicable  to  such  mixtures  at  least  for  data  obtained  on 
charges  with  d  <  7.6  cm. 

The  infinite  diameter  detonation  velocity  for  TNETB  is  given  by 

(mm/psec)  =  1.947  +  3.660  pQ 

and  the  diameter  effect  on  D  for  this  H.E.  is  very  small  in  the  range 
2.54  <  d  <  7.62  cm. 
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APPENDIX  A 

Effect  of  Discarding  Two  Data  Points  of  Table  1 

It  was  mentioned  in  the  text  that  two  values  of  D  at  d  =  6.35  cm 
seemed  out-of-line  in  Figure  1.  When  these  two  points  are  omitted 
from  the  least  squares  treatments,  the  results  of  Table  A1  are  ob¬ 
tained.  The  omission  of  these  points  reduced  the  q.m.e.  of  the  D  vs 
d**1  curve  for  pQ  =  0.722  g/cc  by  a  factor  of  about  3;  it  did  not 
change  the  q.m.e.  of  the  curve  for  pQ  =  0.701  g/cc.  The  q.m.e.  of 
the  resulting  D^  vs  curves  was  reduced  by  a  factor  of  about  2.6. 
The  improvement  is  evident  in  the  resultant  curves.  Figures  A1  (pQ  = 
0.722)  and  A2,  but  does  not  seem  large  enough,  compared  to  the  ex¬ 
perimental  scatter,  to  justify  discarding  the  points  at  d  =  6.35  cm. 
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TABLE  A1 

RESULTS  WITH  OMISSION  OF  DATA  AT  d  =  6.35  cm. 


Pq 

0.722 

0.701 


D^  (mm/ixsec) 


Min. 

2.956  (O.OI67) 
2.910  (0.044) 


Max. 

2.970  (0.0134) 

2.944  (0.041) 


-b 

cm  mm/u.sec 
Min.  Max. 

3.661  3.768 

2.375  2.480 


Minimum 

Maximum 


D±  =  1.057  +  2.620  pQ  (0.0146) 
D±  =  1.100  +  2.505  pQ  (0.0154) 
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APPENDIX  B 

APPARENT  AGEING  EFFECT  OF  AP 

A  small  sample  of  XP-17  was  stored  for  a  year  after  the  original 
study1  was  completed.  The  results  on  charges  (pQ  =  O.685  g/cc)  pre¬ 
pared  from  this  aged  sample  are  as  follows: 


D  mm/nsec 


Shot  No. 

d(in. ) 

Aged 

Original* 

400 

1.0 

*(1.59) 

1.88 

595 

1  5/8 

1.955 

2.10 

592 

2.0 

2.140 

2.28 

*  From  smoothed  curves  of  Ref.  (1), 


X'P-il  after  some  months  storage  became  mildly  caked  and  had  to  be 
rolled  before  it  could  be  used  to  make  charges.  After  caking,  a 
2  in.  diam  charge  at  pQ  =  0.722  g/cc  showed  D  of  2.24  mm/Usec  com¬ 
pared  to  the  precaking  value  of  2.55  mm/|i.sec  (Ref.  2,  smoothed  data). 

Lot  N-128  was  fired  at  pQ  =  0.701  g/cc  and  d  =  2.0  in.  about  three 
weeks  after  it  was  received  and  again  a  month  later.  Results  were: 

Shot  No.  D(mm/u.sec)  Date 

412  2.665  5/51/67 

414  2.497  6/27/67 

None  of  the  changes  noted  are  large  enough  to  be  conclusive  except 
the  failure  of  the  one  inch  diam  charge  of  aged  XP-17.  However,  all 
changes  are  in  the  direction  of  lowered  D  with  Increased  age. 

Future  work  will  explore  this  trend,  particularly  at  low  densities 
and  with  freshly  ground  material,  conditions  which  might  be  expected 
to  favor  the  largest  charges. 
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APPENDIX  C 

DATA  PROM  PRELIMINARY  WORK  ON  AP/HMX  MIXTURES 

Since  supplies  of  HMX  of  various  particle  size  were  on  hand,  explora¬ 
tory  work  on  AP/H.E.  mixes  was  carried  out  with  this  H.E.  On  the 
basis  of  this  work,  as  noted  in  the  text,  it  was  found  that  more 
uniform  charges  could  be  prepared  from  the  80/20  mixes  than  from  the 
90/10  although  both  showed  some  agglomerates  of  the  fine  H.E. 

Table  Cl  contains  all  data  obtained  in  the  preliminary  work.  Com¬ 
parison  of  shots  471  and  472  with  475  and  474  indicates  that  a  fine 
H.E.  is  more  effective  than  a  coarse  in  extending  tne  detonability 
range  of  the  AP/HMX,  90/10,  mixtures;  hence  fine  H.E.s  were  used  for 
the  remainder  of  the  work.  However,  it  should  be  noted  that  (1)  the 
coarser  HMX  mix  could  be  compacted  to  a  slightly  higher  #  TMD  (l# 
higher)  them  the  finer  HMX  mix,  an  effect  of  coarse  particles  fre¬ 
quently  observed  in  compaction,  and  (2)  the  higher  %  TMD  rather  than 
the  coarser  HMX  might  be  responsible  for  the  failure  of  the  2.54  cm 
diameter  charges  to  detonate  at  pQ  *  I.765  g/cc  (90.8^  TMD). 

Where  the  raw  data  can  be  compared,  the  finer  AP  resulted  in  a 
slightly  higher  D  at  the  same  charge  diameter  and  a  slightly  lower 
density;  compare  shot  458  with  464  and  457  with  462.  For  this 
reason  as  well  as  in  hope  of  obtaining  better  dry  mixing,  fine  AP 
was  used  in  the  subsequent  work. 

The  overall  gradient  of  the  D  vs  p_  curve  for  2.54  cm  diameter  charge 
of  AP/HMX,  90/10,  is  about  1.64  mm  p,sec  cc  g  .  (See  Figure  Cl). 
This  is  about  comparable  to  the  slopes  found  for  AP  alone  (1.8  to 
2.28  for  d  of  5.81  to  7.62  cm)  As  pointed  out  above,  however,  the 
high  density  point  may  be  close  to  the  failure  limit,  and  consequently 
the  curve  may  not  be  linear  over  the  entire  range.  It  should  also  be 
noted  that  this  AP  alone  is  not  detonable  in  a  2.54  diameter  charge 
at  any  density,  and  we  have  been  able  to  detonate  our  rinest  AP  only 
up  to  about  81#  TMD. 

Data  for  AP,N126,  which  is  slightly  finer  than  AP,N127. 
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APPENDIX  C  (Cont'd) 

Non-linearity  is  exhibited  by  the  D  vs  pQ  data  of  the  80/20  mixture 
(d  =  2.54  cm)  also  plotted  in  Figure  Cl.  The  slope  is  approximately 
4.43  in  the  range  87.5-92.4#  TMD  but  the  curve  appears  to  bend  over 
between  92.4  and  95.9#  TMD.  This  might  be  experimental  error  in  the 
last  point,  or  a  high  density  behavior  comparable  to  that  of  TMT, 
high  bulk  density  NQ,  or  AP  itself.  In  any  event  addition  of  HMX  to 
AP  extends  its  detonability  in  2.54  cm  diameter  charges  from  about 
51#  TMD  to  at  least  88.8#  (10#  HMX)  and  95.9#  (20#  HMX). 

There  are  also  data  for  three  D  vs  d  series  in  Table  Cl.  These  are 
plotted  in  Figure  C2  where  each  of  the  three  curves  shows  the  break 
indicating  an  oxidation-reduction  reaction,  as  discussed  in  the  text. 
Although  the  finer  AP  appears  to  produce  larger  D  at  any  given  pQ  and 
d,  it  doe3  not  seem  to  affect  the  value. 

The  additivity  rule  was  not  applied  to  these  results  because  it  is 

not  quantitatively  descriptive  of  porous  mixtures.  The  additivity 
24 

rule  gives  the  ideal  detonation  velocity  of  a  mixture  as 


°iv  x3  \-}  <«> 

where  x  is  mass  fraction  and  the  subscripts  i,  j  and  v  refer  to  ideal 
(Infinite  diameter),  component  j,  and  voidless  density,  respectively. 
At  100#  TMD,  Eq.  (Cl)  gives  the  detonation  velocity  of  cyclotol 
mixtures  to  0.3#.  When  an  analogous  approximation  is  used  for  porous 
charges,  the  component  and  mixture  velocities  must  be  computed  at 
the  same  #  TMD,  but  the  approximation  becomes  poorer  as  the  #  TMD 
decreases.  Thus  at  70 #  TMD,  the  Di  computed  for  the  mixture  from 
known  ideal  velocities  of  the  components  is  about  1#  low.  The  error 
increases  at  lower  #  TMD, 

When  this  approximation  is  used  to  compute  the  velocity  of  one  com¬ 
ponent  from  the  measured  velocity  of  the  mixture  and  known  velocities 
for  other  components,  the  error  is  roughly  twice  that  of  the  inverse 
procedure.  Hence  our  effort  to  obtain  mixtures  of  #  TMD  >  95  before 
applying  the  additivity  rule. 
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1.7 


1.9 


DENSITY  p0  (g/cc) 

FIG.  Cl  D  VS  p0  CURVES  FOR  AP/tHMX  AT  d  =  2.54  cm 
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PRELIMINARY  UNCORRECTED  SHOT  DATA  ON  MIXTURES  OP  AP/HMX 
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a  Corrected  to  pD  =  1.576  g/cc 
b  Corrected  to  pQ  -  1.576,  same  elope  used. 
•*.  Corrected  to  p0  «  1.74  g/cc 
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